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Abstract: We propose a new doorway to study the interplay between equations of state 
of dense matter and compact stars in gauge/gravity correspondence. For this we con- 
struct a bulk geometry near the boundary of five-dimensional spacetime. By solving a 
constraint equation derived from the bulk equation of motion together with the Tolman- 
Oppenheimer-Volkoff equation, we determine the equations of state for compact stars. The 
input parameters in this study are the energy density and pressure at the center of the 
compact objects. We also study how the equation of state depends on the parameters. 
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1. Introduction 

The equation of state (EOS) of dense matter finds its pivotal importance in various physical 
systems such as compact astrophysical objects and atomic nuclei. Using a given EOS of 
dense (hadronic) matter, we can determine the mass and radius of a neutron star by 
solving the Tolman-Oppenheimer-Volkov (TOV) equations. Or, using the observed mass 
and radius of a neutron star, one can constrain the EOS of dense matter ||, ||] . Ultimately, 
one can explore the one-to-one correspondence between an EOS of dense matter and a 
mass-radius relation of a compact star. 

Recently there have been several interesting attempts based on AdS/CFT j|, ||, ||| 
to understand compact stars ||, pXj] , self-bound dense objects such as atomic nuclei 

12, [l4|], and a two-dimensional plasma ball fl5| . 
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In this work, we try to understand the compact stars in gauge/gravity duality. Since 
one essential feature associated with the compact star is the gravitational force, one has to 
construct field theory coupled to gravity. However, the AdS/CFT is achieved by taking the 
decoupling limit where the gravity is completely decoupled from a boundary field theory 
system. Therefore, one has to properly introduce the boundary gravity to any constructions 
for gravitationally bound compact stars. Such a construction can be realized through an 
effective action 

Seff = Wcft{9{o)] + 16 l G J <? x y/-9(p) R ho)] + ■■■ , (1-1) 

where Wcft [<7(o)] i s the generating functional of the connected Green's functions for con- 
formal field theory (CFT) with the boundary metric g$ and G4 is the four-dimensional 
Newton constant. The dots denote higher derivative curvature terms. A well-established 
way to obtain an effective action of a boundary field theory system coupled to gravity in 
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AdS / CFT is to introduce a finite UV-cutoff to the CFT system 1 . The cutoff defines an effec- 
tive four-dimensional Newton constant along with other parameters from five-dimensions. 



This generalization of the AdS/CFT correspondence has been studied in |17], 18], see [19] 
for a recent review. We will provide this kind of example in appendix |A|. 

With the generalized AdS/CFT (or gauge/gravity) duality at hand, we try to build up 
a holographic model in which conformal field theory interacts with the boundary gravity 
and possibly uphold a stellar object. For this, we take a simple ansatz for the bulk metric 
and the bulk matter field. Then, we solve the bulk equations of motion near the boundary 
in Fefferman-Graham coordinates. Our strategy to deal with the holographic perfect fluid 
stars is to put perfect fluid star metric as a boundary metric and find the corresponding 
bulk geometry near the boundary of five-dimensional spacetime. In this procedure, we 
obtain a constraint from the bulk equation of motion. The constraint together with the 
TOV equations will determine the EOS of a dense matter. As it is, this observation 
stands quite opposite to a standard method in which an EOS is an input to solve the 
TOV equations. Mathematically, it is quite obvious that we could fix the EOS with the 
additional constraint from the bulk. However, physically it is not clear how the constraint 
from the extra dimension could affect the EOS in four-dimensional spacetime. With our 
results at hand, we study the interplay between equations of state of dense matter and 
perfect-fluid stars and classify sets of our model parameters. The values of our model 
parameters are specified at the center of a the perfect-fluid star. We also consider a case in 
which we impose boundary conditions at the center and at the surface of the perfect-fluid 
star. 

This paper is organized as follows. In section ^, we briefly review compact stars. In 
section |3|, we consider the simplest case in a holographic approach, the pure gravity system 
in the bulk, and we explain our assumptions and ansatz for the perfect-fluid star. In section 
we generalize the simplest case by adding a bulk scalar field. We then show that this 
modification makes more realistic configuration possible. In section ||, we make a brief 
summary of our present study and provide some future extensions. In appendix 0, we 
review a construction for conformal field theory coupled to gravity studied in 



2. Compact stars in a nutshell 

The formation of a (gravitationally bound) star is mainly governed by two forces: gravita- 
tion and thermal/degeneracy pressure generated by the matter inside the star. In the case 
of compact stars, the degeneracy pressure is dominated by neutrons (neutron stars) and by 
electrons (white dwarfs). Note that the energy of the white dwarfs is dominated by nuclei. 
The force balancing between the gravitational collapse and the internal pressure leads to 
a basic structure equation for stars, see figure [I], 

AdP = - G4m( r r 2 )dm , (2.1) 

1 There is an alternative way suggested recently, where normalizable gravity wave functions are introduced 
by an IR-cutoff |16|. 
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Figure 1: Sketch for hydrostatic equilibrium. Here, F denotes the gravitational force and AdP is 
for the force from the thermal/degeneracy pressure. 



where AdP denotes a net force on the mass element due to the pressure difference dP. 
Since dm = p(r)Adr, we obtain the equation of hydrostatic equilibrium, for instance see 

dP(r) G^m(r)p(r) 



dr r 2 

i2. 



(2.2) 



m(r) is the mass interior to a radius r, m(r) = J*J* p(r')4mr' dr'. From this, we derive 

dm(r) 



dr 



Airr 2 p{r) . (2.3) 



Note that we do not distinguish the matter density from the energy density (= matter 
density x c 2 ) since we take c = 1, where c is the speed of light. 

The treatment in the above is valid in the Newtonian limit. A gravitational effect, 
however, is well-described by general relativity. So we move onto a solution of the Einstein 
equation that is essential to study a static star with spherical symmetry, see [21| for more 
details. For a static spherically symmetric stellar object, the natural and obvious metric 
ansatz is 

ds 2 = _ e mr) dt 2 + e 2h(r) dr 2 + r 2 dQ 2 ^ ^ 

where dfi 2 . = dO 2 + s'm 2 6d<p 2 . With this ansatz we solve the four-dimensional Einstein 
equation 

R/xv ~ ^9fivR = 87rG 4 Tfa, . (2.5) 

Here, we make an assumption that the matter inside the star is perfect fluid. Then, the 
energy-momentum tensor can be written as 

T^ = (p + P)/u" + Ff , (2.6) 

where is the 4- velocity of the fluid with the normalization: u^u^ = —1. From the 
(t, t)-component of the Einstein equation, we obtain 



1 d 

r 2 dr . 



r (1 -e- 2/l(r) )] =8nG 4 p{r). (2.7) 
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With 2Gim(r) = r(l — e 2/l ( r )), we can rewrite this equation as 

4Trr 2 p(r), (2.8) 



dm(r) __ t _. 



dr 

which is the same with (|2.3l). Then, the metric (|2.4| ) changes into 



ds 



2 



-e^dt* + (l - ?^T\ 1 dr 2 + r 2 dQ 2 (2 _ 



On the other hand, from the (r, r)-component of the Einstein equation, we obtain 

G 4 (m(r) + 47rr 3 P(r)) 
= V- 2G 4 m(r)) ■ (2 ' 10) 

The relativistic extension of hydrostatic equilibrium is derived from the conservation law 
of energy-momentum tensor. 

V M T^ = => P'(r) = -(P(r)+p(r))f'(r) (2.11) 

Inserting (2.10) into above equation, we arrive at the TOV equation 

dP{r) _ G A (m(r) + 4vrr 3 P(r)) (P(r) + p{r)) 
dr r (r — 2G4m(r)) 

which is the general relativistic extension of (|2.2| ). To solve the equations, (|2.2[) (or ( |2.12|) ) 
and (^.3|), one generally needs to impose two initial conditions at r = 0, m(0) = and 
p(0) = p c , and the EOS, P = P{p). Usually one has to solve the equations numerically, but 
there are several analytic solutions. The simplest analytic solution is the uniform density 
star, p(r) = p c . In this case, the corresponding inner pressure is given by, see [21] for 
instance, 

w^^-^/^-d-^/^ (213) 

3(1 - 2G 4 M/R) 1/2 - (1 - 2G 4 Mr 2 /R 3 ) 1/2 
Another interesting analytic solution is so called Tolman VII solution: 

p(r) = p c (1 - (r/Rf) , (2.14) 

where p c is the central density and R is the radius of the star. The corresponding pressure 
is given as an analytic form. Since it is too complicated, we will not show the pressure 
here. For a summary of the analytic solutions, we refer to HJ. 

Outside of the star, the density p{r) and pressure P(r) vanish since there is no matter. 
Then, from the equation ( [2 .8; ) we can see that m{r) becomes a constant. When we integrate 
the equation fl2.10| ) with the condition P(r) = and m(r) = m(R) = M, it requires a 
boundary condition /(oo) = since the metric should be asymptotically flat. Then, the 
final form of the exterior metric is 



ds 2 = -ll — Jdt 2 +n — J dr 2 +r 2 dQ, 2 2 . (2.15) 

This is known as the Schwarzschild metric. From now on, we will set G4 = 1. 
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3. Uniform density stars 



We start by considering the five-dimensional vacuum Einstein equation with a negative 
cosmological constant. 

1 6 

Rmn ~zGmnR "j^Gmn = (3-1) 

From the result in p2| , one may bring the general behavior of the bulk metric in the 
Feffermann-Graham coordinate as 

L 2 

ds 2 = g mn dx m dx n = (dz 2 + g^(x, z)dx fl dx v ) , (3.2) 

where 

CO 

g^x, z) = $(x) + z 2 g$(x) + z*g$(x) + £ *g${x) . (3.3) 

n=5 

We confirmed that the coefficients of z and z 3 are identically zero. In general, there can 
be a logarithmic term in this expansion. This logarithmic term contributes to a scheme 
dependent part of the boundary energy-momentum tensor, and we will neglect this term 
for simplicity. Mathematically, this choice restricts g^ix, z) to be the polynomial function 
of z near the boundary of the AdS space. Since we are going to describe a static spherically 
symmetric stellar configuration, we take the following ansatz for the metric. 

g^(x, z)dx fl dx u = g tt (r, z)dt 2 + g rr (r, z)dr 2 + g gg (r, z) (d9 2 + sin 2 Odcj?) , (3.4) 

where r is the radial coordinate of the boundary space. 



Now, we are ready to solve the bulk Einstein equation (3.1) near the boundary. Putting 

(n) 

the metric ansatz into the Einstein equation, we obtain differential equations of g^J with 
respect to r. In the leading order of the expansion near the boundary, O(z ), all the 
non-trivial components of gfy can be determined in terms of gffi 



This is the same result with the one in p2| , p3[ . The equations obtained in the next- 
to-leading order, 0(z 1 ), are trivially satisfied with the relation in ( |3.5D . Since we are 
interested in a starlike object of the boundary system, we need to deform the boundary 
metric. Without loss of generality, we consider the static spherically symmetric metric 
discussed in the previous section, 

9i? = , ^ = (l - ^) " and g$» = r 2 . (3.6) 

Using this ansatz, the expression in ([T5|) is given as follows. 

9tt ( r ) = — {m'{r) - (3m(r) + r (-2 + m'(r))) f'{r) 

+r (r - 2m(r)) f{r) 2 + r(r - 2ra(r))/"(r)) (3.7) 
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giV(r) = — tt, ttt (3m(r) — 2rm!(r) — r (r — 3m(r) + rm'(r)) /'(r) 

3H (r — 2m(r)) / 



+r 2 (r - 2m(r)) (/'(r) 2 + f"{r))) (3.8) 

( .SmlV") -I- rmf(r) — r (r — 3?i 

6r 



9eei r ) = (3m(r) + rm'(r) — r (r — 3m(r) + r?n'(r)) /'(r) 



+r 2 ( r _ 2m ( r )) (/'( r )2 + /"( r ))) (3.9) 

Considering the equations in the next-to-next-to-leading order, 0(z 2 ), we obtain a con- 
straint for the boundary metric. The constraint has a little bit complicated form: 

= — 3m(r) 2 + r (— 4r + 6m(r) + rm'(r)) m'(r) + 2r 2 (r — 2m(r)) m"(r) — 2r (2r 2 

+ (— 5r + 3m(r)) m(r) + (— 3r + 4m(r) + rm'(r)) m'(r) + 2r 2 (r — 2m(r)) m"(r)) /'(r) 

+ r 2 (9r 2 + (— 32r + 29m(r)) m(r) + r (— 4r + 6m(r) + rm'(r)) m'(r) 

+2r 2 (r - 2m(r)) m"(r)) /'(r) 2 + 2r 3 (r - 2m(r)) (-3r + 5m(r) + rm'(r)) /'(r) 3 

+ r 4 (r — 2m(r)) 2 f'{r) A — (4r 2 (r — 2m(r)) (— r + m(r) + rm'(r)) 

-4r 3 (r - 2m(r)) (-m(r) + rm'(r)) /'(r) + 2r 4 (r - 2m(r)) 2 /'(r) 2 ) /"(r) 

+ r 4 (r - 2m(r)) 2 f"{rf - 2r 3 (r - 2m(r)) 2 (-1 + rf'{r)) f"(r) . (3.10) 

In general, it is not easy to solve this equation since it is hard to impose suitable boundary 
conditions. To make this equation in a more tractable form, we constrain our system 
further to consider perfect fluid stars. In order to describe them, we put the following 
restriction on our metric as briefly summarized in the previous section, see [21] for details. 

m{r) + 47rr 3 P(r) 

f(n = — / 7 vi » 3 - n 

r (r — 2m{r)) 
rr 

m(r) = 4-7T / p(r')r' dr' , (3-12) 



„// s . . . . NN m(r) + 47rr 3 P(r) 

P'{r) = - P r + P r 4 9 ^ . 3.13 

where /o(r) and P(r) are the energy density and pressure of a perfect fluid star. Exploiting 
the above equations, we rewrite ( 3. 1C| ) in a much simpler form 

= (-P(r) + (3m(r) - 4irr 3 p(r)) 
' ' r {r — 3m(r) — 4-7rr 3 P(r)) 

Now one can numerically solve ( 3.12j ), ( |3 . 1 3| ) and ( p.!4| ) with the following boundary con- 
ditions: the value of the pressure and the energy density at the center of the perfect fluid 
star. Note that m(0) = by construction. An example of such solutions is given in figure |2|. 
The example is a well-known uniform density star of the four-dimensional Einstein gravity 
in flat spacetime. This configuration does not show interesting structure. However, it is 
not a meaningless result in the following sense. This simple analysis demonstrates that our 
five-dimensional bulk construction near the boundary does constrain the four-dimensional 
gravity system, making it possible to determine the EOS and to study properties of the 
perfect fluid star. 
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Figure 2: Energy density and pressure of a uniform density star 

Before we move to a realistic case, we discuss a bit more on the uniform density 
solution. Though it is difficult (or almost impossible) to solve the coupled non-linear 
differential equations analytically, we can study the behavior of the uniform density solution 
analytically near r = thanks to the regularity condition at the center. From the relation 
( |3.12 ), p(r) can be replaced with m{r) as p(r) = m' (r) / (Atrr 2 ) . Furthermore, considering 
the regularity of p(r) at r = 0, one can impose the boundary conditions, m'(0) = m"(0) = 0. 
Now, we will expand P(r) and p(r) through Taylor series around r = 0. Here, we focus on 



= 



r r 



- 3m(r) 

- 3m(r) 
m( 4 )(0) 2 



r [r 



4irr 3 P(r)) p'(r) - (P(r) + p(r)) (3m(r) - 4irr 3 p(r)) 

s\ d ( m'(r)\ /„ . . m'(r)\ , , . 
4rr P ^ Tr ( S) ~ ( P(r) + 4^ J (MO - 



rm 



(5) 



(0) r 3 



(40vrP(0) + 5m( 3 )(0)) m,( 4 )(0) - 2m( 6 )(0) 



+ ... (3.15) 



24vr ' 48vr 320vr 
To satisfy the equation, it is clear that m^(0) = m^(0) = 0, also m^(0) = 0. Inserting 
these values into the original expression and performing the similar process repeatedly, we 
can see the higher-order derivative values of m(r) at r = vanish except m^ 3 \0) = 87rp(0). 
This almost excludes possibility that all derivative terms of p(r) do not vanish at r = 0, 
and thus we may take p(r) = p(0) = const. Here the word "almost" means that we didn't 
check if all the n-th higher-order derivative of m{r) vanish at r = 0. 

4. Non-uniform density stars 

Since the previous result is too simple to describe realistic compact stars, we try to improve 
the situation by introducing a bulk matter field. It is natural to expect that this matter 
field will lead non-uniform energy density since it will contribute to a r-dependent boundary 
action. It is obvious from ( 3.13| ) that the pressure will change (decrease) with the radius 
r. For simplicity, we consider a bulk scalar 2 and choose its bulk mass as m 2 = — 3/L 2 . 
The Einstein equation we have to solve is now 



R„ 



1 



QmnR 



6 



(4.1) 



2 In a usual holographic QCD, this bulk scalar field is dual to a quark bilinear operator. 
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1 ~ , „ , 1 „ / . .« -2 . 2-2 



r mn = -9 m </»a„</» - -g mn \jd<t>y + m%4> 1 ) • (4.2) 

Here we put the five-dimensional gravitation constant G§ into the scalar field eft. In addition, 
we should solve the following equation of motion for the bulk scalar field 

(V 2 -m 2 )<A = 0. (4.3) 

Now, we take an ansatz for the bulk scalar field. 



[r,z) = z _ 

n=3 



1 (r) + J2z n Mr) (4.4) 



Similar to the metric expansion, the logarithmic term is neglected. The form of the metric 
is the same with ( |3.3| ) near the boundary. Using these ansatz, one can solve the equations, 
Q4.1| ) and ( |4.3| ), perturbatively for small z. For the matter part, we obtain ( [4.5| ) for <f>\ and 
flOD instead of C ^TTOj) . 



= r 2 ^i(r) 3 + (2m' (r) + (-2r + 3m(r) + rm'(r) - r (r - 2m(r)) /'(r)) /'(r) 

— r (r — 2m(r)) /"(?")) <^>i(r) + 3 (— 2r + 3m(r) + rm (r) — r (r — 2m(r)) f'(r)) <j>'\{r) 

- 3r (r - 2m(r)) 0"(r) (4.5) 

= i r 6 (/>i(r) 4 + 2r 3 (m(r) - r (r - 2m(r)) /'(r)) 0i(r) 2 

- 2r 4 (r - 2m(r)) (2 + rf'(r)) <pi(r)<f>[(r) - 3r 5 (r - 2m(r)) ^i(r) 2 

- 3m(r) 2 + r (— 4r + 6m(r) + rm'(r)) m'(r) + 2r 2 (r — 2m(r)) m"(r) — 2r (2r 2 

+ (— 5r + 3m(r)) m(r) + (— 3r + 4m(r) + rm'(r)) m'(r) + 2r 2 (r — 2m(r)) m"(r)) /'(r) 

+ r 2 (9r 2 + (— 32r + 29m(r)) m(r) + r (— 4r + 6m(r) + rm'(r)) m'(r) 

+2r 2 (r - 2m(r)) m"(r)) f'(r) 2 + 2r 3 (r - 2m(r)) (-3r + 5m(r) + rm'(r)) f'(rf 

+ r 4 (r — 2m(r)) 2 /'(?*) 4 — (4r 2 (r — 2m(r)) (— r + m(r) + rm'(r)) 

— 4r 3 (r — 2m(r)) (— m(r) + rm'(r)) /'(?*) + 2r 4 (r — 2m(r)) 2 f'(r) 2 ^j f"(r) 

+ r 4 (r - 2m(r)) 2 /"(r) 2 - 2r 3 (r - 2m(r)) 2 (-1 + rf'(r)) f"'(r) (4.6) 

As we did in the previous section, we consider perfect fluid star configurations, and so we 
adopt the relations in ( 3.1l| ), (|3.12|) and (|3.13| ). Then, the two equations, ([4.5|) and ([O]), 
become much simpler, 

^ = 5 i \ i ^ Mr) 3 ~ 47rr 2 0!(r)(3P(r) - p{r)) 
Sr (r — 2m[r)) 

- 6 (r - m(r) + 2vrr 3 (P(r) - p(r))) <^(r)) , (4.7) 

+16nrP(r)4> l (r) [<f>i(r) + n^r))) + 32vr (-3m(r) + 4vrr 3 P(r)) p(r) 
+12m(r) (-8vrP(r) + r^i(r) (^i(r) + r4>\(r))) + 128vr 2 r 3 p(r) 2 ) . (4.8) 
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To solve the equations, (|3.12[) , (pL"7[) , and ( @g ), we need to impose three boundary 

conditions at r = 0: p(0) = pc, -P(O) = P c and <^>i(0), 3 together with regularity condition at 
the center of the object. Here 4>i(0) has no significant physical meaning at this stage, but 
it may be associated with the energy density on the surface of our spherical objects, see 
below for more on this. An example of the solutions is given in figure ^, where we choose 
p c = 2.5713 x 10 9 MeV 4 , P c = 3.2141 x 10 8 MeV 4 and ^i(O) = 0.2703. In order to have some 
idea about corresponding baryon number density at the center of the star, we resort to the 



pressure-baryon density plot in [24] and find that our chosen P c corresponds to roughly 3.5 



times the normal nuclear matter density. In [24|, based on the recently determined mass- 
radius relation of neutron stars, the pressure of cold dense matter is measured. In this way, 




Figure 3: Energy density, pressure and equation of state for p c — 2.5713 x 10 9 MeV 4 and P c = 
3.2141 x 10 8 MeV 4 with surface energy density 3.26 x 10 s MeV 4 . This configuration provides a star 
of the radius 11.45 km and the mass 1.26 Mq. 



we can obtain various EOSs once we specify the three values at the center of the compact 
stellar object. However, among the three, the value of </>i(0) is entirely arbitrary. To remove 
this arbitrariness and to reduce the number of input parameters, we assume that the energy 
density at the surface p(R) is zero. In this case, we can read off the value of <^>i(0) from 
our solutions. In general, the surface density is much smaller than the central density, and 
therefore our assumption is a plausible leading order approximation. Then, our parameter 
space becomes two dimensional parameterized by the central values of the energy density 
and the pressure (p c ,P c ). Though the physical meaning of our inputs is now clear, we 
still have an ambiguity since we have no rigorous way to relate p c to P c . For illustrative 
purposes, we choose P c = ap c with a = 1/3, 1/4, 1/5, 1/6, 1/8. The resulting mass-radius 
relation is summarized in table |]. For every point in the parameter space, we can find 
the corresponding equation of state. This means that a set of central density and pressure 
defines a unique equation of state. For instance, we demonstrate the energy density, the 
pressure and the equation of state for p c = 2.5713 x 10 9 MeV 4 and P c = 3.2141 x 10 8 MeV 4 
in figure In this case, we obtain </>i(0) = 0.2704 from our solutions. 

Now we try to understand the physical roll of </>i(0). If we change the value of </>i(0) 
above, we will obtain a star configuration with different surface energy densities. Thus, 
with given values of the central density and pressure as the boundary condition, (f)\ (0) is 
closely tied to the surface energy density. This does not necessarily imply that <^i(0) is 



Since (4.7) is the second order differential equation, we need to have one more boundary condition, 



4>'i(0)- However, we can express it in terms of the other three conditions due to the regularity condition at 
the center. 
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Lvp 


Mass (M ) 


1.59 


1.38 


1.23 


1.12 


1.04 


0.97 


0.92 


0.87 


Radius (km) 


12.33 


10.68 


9.55 


8.72 


8.07 


7.55 


7.12 


6.75 


Pc = gPc 


Pc 


op 


Ap 


hp 


6p 


7p 


Rp 


9p 


Wp 


Mass (M ) 


1.26 


1.09 


0.97 


0.89 


0.82 


0.77 


0.73 


0.69 


Radius (km) 


11.66 


10.10 


9.03 


8.25 


7.64 


7.14 


6.73 


6.39 



Table 1: Sample mass-radius relations. Here p = 8.5704 x 10 8 MeV 4 . 



physically tied to the surface energy density. For instance, if we use the surface energy 
density and the central pressure as the boundary condition, then 4>i(0) will be related 
to the central density. Therefore, within our framework, it is hard to put some physical 
meaning on </>i(0). In terms of the holographic QCD, one may say that <?H(0) is nothing 
but the quark mass at the center of the star by the AdS/CFT dictionary and 4>\{r) is 
an effective density-dependent quark mass since the density varies with the radius of the 
compact star. 

Before closing this section, we give a remark on the EOS obtained in this study. In 
a usual approach to calculate the mass-radius relation, one adopts an EOS from various 
studies in flat spacetime with no general relativistic effects and then solve the TOV equa- 
tion. Though it is obvious, we remark here that this standard approach, TOV equation 
with flat EOS, is perfectly consistent with general relativity. In our case, we solved the 
TOV equation together with the constraint from the bulk to obtain the EOS. In this sense 
we can argue that our EOS includes the effect of general relativity (or gravity). 




4 6 8 10 2 4 6 ft 10 5.0x10 s l.OxlO 9 UxlO 9 2.0xl0 9 2.5 x 10 9 



Figure 4: Energy density, pressure and equations of state for p c = 2.5713 x 10 9 MeV 4 and P c = 
3.2141 x 10 8 MeV 4 . This configuration gives the radius 11.66km and the mass 1.26 M@. 
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5. Summary 

In this work, we investigated stellar configurations in gauge/gravity correspondence to 
study perfect fluid stars. For this, we have considered holographic dual of a boundary field 
theory system coupled to gravity. To retain the gravity degrees of freedom, we introduced a 
finite UV-cutoff that defines the four-dimensional Newton constant. With the assumption 
that the energy density at the surface is zero, we determined the EOS and mass-radius 
relation of a compact star, once the central values of the energy density and pressure are 
specified. An ambiguity in this case is that the central values of the energy density and 
pressure are not correlated each other. In this study, the EOS is not an input to solve the 
TOV equation, but the output of our analysis together with the corresponding mass and 
radius of a star. 

So far, we analyzed the bulk metric only near the boundary. To explore the full 
bulk metric, we have to solve partial differential equations with appropriate boundary 
conditions. Our study could be very important information for the boundary data. In 
future, we will seek a complete solution for the bulk metric that covers whole bulk space. 
Another possible extension is to generalize our ansatz. Since the generalization can affect 
directly energy-momentum tensor, stars can have interesting structure even in the pure 
bulk gravity case. A study with a more general ansatz for the bulk metric is relegated to 
a forthcoming publication. 
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A. AdS/CFT with gravity 

In this section, we provide an explicit example of a generalization of the AdS/CFT cor- 
respondence with gravity. In []l7| , |lg[| , it is shown that the well-known RS-II model pl| 
is dual to a system which is similar to ( |1.1| ). Here we present a compact review on the 
derivation in jlSj] . 

The AdS/CFT correspondence is defined by equivalence between the partition function 
of five-dimensional gravity theory and the partition function of four-dimensional conformal 
field theory as follows. 




(A.l) 
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where 7 is the induced metric at the boundary of AdS space and WcftM denotes the 
generating functional for connected Green's functions. 5Q and ip are path-integral measures. 
The action for the five-dimensional gravity theory is given by 



' gravity 



16vrG 



l — [ d 5 xV^Q(R + 2A) + — L- / (fx^ K + 5i + S 2 + S 
tG 5 J 87tG 5 J 



where the five-dimensional metric can be decomposed into 



ds 2 = g mn dx m dx n = G^(x, z)dx^dx u + 



L 2 dz 2 



(A.2) 



(A.3) 



and the range of coordinate z is zq < z < 00 with zq <C L. The UV-cutoff A is related 
to zq (~ •jj). The first term of ( A.2Q is the the well-known Einstein-Hilbert action with 
the cosmological constant, and the second term is the Gibbons-Hawking term. The other 
terms are introduced by holographic renormalization as follows. 4 



Si 
S 2 
S3 



8vrG 5 L 
L 



d 4 x, 
d 4 x x 



32vrG 5 
L 3 log(z /L) 
64vrG 5 



-7, 



R^(9(0))R^(9(0))-^R(9(0)f 



(A.4) 
(A.5) 
(A.6) 



These are counter terms needed to cancel divergence in the Einstein-Hilbert action and 
the Gibbons-Hawking term. On the other hand, in view of four-dimensional boundary 
theory, these terms can be interpreted as follows. Si is just four-dimensional cosmological 
constant, and S2 is four-dimensional Einstein-Hilbert action. S3 could be interpreted as 
higher derivative corrections with coupling that depends on the cutoff scale. Actually, it is 
proportional to the conformal anomaly. 

Now, let us turn to the RS-II model, which is defined as, 



Srs-ii = Seh + Sqh + 2Si + S n 



(A.7) 



where Seh an d Sqh are the Einstein-Hilbert action with cosmological constant and the 
Gibbons-Hawking term as in the previous case. In addition, a brane exists at z = zq and 
the 25i gives its tension, j • The model has confined matter on the brane at z — z§ 1 
which is described by the action S m . Let's consider path-integral of this action. One can 
write down the partition function as follows. 



Zrs-ii[i] 



iSrs-ii [Go+SQ, Sip] 



d[5Q]d[5^\e 



2iSi[7] / d[SQ]d[5ip]e iSEH ^ 0+s ^ +iSGH ^ 0+5 ^ +iSm ^' fi '' y ^ 



(A.8) 



3 2tSi[ 7 ] 



d[5G]e 



iS EH [Go+SG] +iS GH [Go+SG] 



one— bulk 



d[8ip]e 



iS m [5tp;y] 



Where gj^J denotes z-independent part of (z 2 / L 2 )Q l _ lv 
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Here to arrive at the final expression, we took account of the Z 2 symmetry, two bulk- 
regions, of the RS-II model. Using the fact that integration of Seh and Sqh already 
appeared in the ( |A.1[ ) 5 , we can express the partition function as follows 

Zrs-iiH = e 2i (^FT-5 2 -5 3 )[7] J d[S(p]e iS ™^^ . (A.9) 

As mentioned earlier, S m is the action for the matter field confined on the brane in the 
RS-II model, thus if we ignore the matter action, the extremum of the RS-II model is 
described by following effective action, 

Seff = WcftH - S 2 [j] - 5 3 [ 7 ] • (A.10) 

Now, it is obvious from this equation that the RS-II model is dual to the conformal field 
theory coupled to gravity including higher curvature interaction. For more clear application 
to our paper, we use the result of J22| |. Solving Einstein equation in five-dimensional AdS 
space, the above induced metric G/j,u(x,z) can be decomposed as follows. 

t2 ( 00 \ 

^ = 72 U$V) + z2 9tK*) + z'\ogz 2 h, v (x) + (A.ll) 

If gffi and gffi are given, the other terms are determined. Here gffi is interpreted as 
the boundary metric at the leading order in ^ and gfy is closely related to boundary 
energy-momentum tensor. Then we may write down above effective action in terms of the 

L 3 



boundary metric gffi as follows. 



Seff = Wgft[9 { o)] + 32 ^ G5Z 2 J *XyFm R \9Vti ( A " 12 ) 



where ". . ." denotes higher curvature terms from S^fy]. 

Let us consider the four-dimensional Einstein equation derived from this effective ac- 
tion. After variation, we can easily expect the following result 



R 



[9(0)} ~ \g { $R\9<f»\ = 8vrG 4 {T^ FT + H^ + ... ) , (A.13) 



where T^ U FT and come from the variations of Wcft and logzo term with respect to 
g$ and ". . ." means the contribution coming from the higher derivative correction terms. 



Thus we may regard the right hand side of ( A.13| ) as an effective energy momentum tensor 



T^J = T£f T + H^ v + . . . , where H^ u is proportional to h fMU (x) in QA.11]) . 



So far, we have reviewed part of [18] and translated their results in our convention. 
From the result, we can read four-dimensional gravitational constant as G4 = 2G^Zq/LP' 
from ( A.12| ). Note that in the limit zq — * 0? G4 vanishes and gravity is decoupled from the 



field theory system. 



S 1 +iS 2 +iS 3 f r l\fiC]p iSEH+iS G H = e iW CFTh) 



5 Z[ 7 ] = e xSl+iS2+lS3 f d[8G] 
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